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1 Introduction 

In the recent past the success of percolation models 
and their link with other generic approaches like Ising and 
Potts 
models 



has led to the development of many lattice 
which were used as sensible albeit schematic 
descriptions of excited disassembling nuclei. As simple as 
they may appear, their thermodynamic properties were 
considered as being at least qualitatively those of bound 
nucleon systems which interact essentially by means of the 
short range nuclear interaction. It was tacitly implied that 
quantum effects do not qualitatively alter those properties, 
at the excitation energies which characterize fragmenting 
nuclei. 

The introduction of time-independent stationary de- 
scriptions of fragmented systems presupposes that these 
systems are in thermodynamic equilibrium. Whether this 
is a realistic assumption and is effectively realised is still 
an open question. There exist however many data which 
show that most experimental results can be understood in 
this framework ||^. 

The present contribution follows a double objective. 
First, in the framework of lattice models, we want to anal- 
yse the effect of the Coulomb interaction which is superim- 
posed onto the short range nearest neighbour interaction 
which mimics the nuclear potential. This is done in the 
framework of the so called Ising Model with Fixed Mag- 
netization (IMFM) Ijlffl and its extension to the case where 
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proton-neutron are differentiated from proton-proton and 
neutron-neutron interactions. The treatment of the long 
range interaction is performed on the same footing as the 
short range one i.e. without any approximation in the cal- 
culations of the interaction between a given proton and all 
the others in the system. In a second part we relax the lat- 
tice structure of the system. Fragmenting nuclei are disor- 
dered systems whose constituents are not located at fixed 
positions on a regular lattice like in crystals but more in 
the continuum of position space. A priori, a more realistic 
description of such systems is realised in the framework of 
so called cellular models . We show and discuss whether 
the freedom to occupy any space position leads or not to 
qualitative changes in the thermodynamic and topological 
(fragment formation) properties of the system. 

In section 2 we present a sketchy description of the 
IMFM model and the way we implement the Coulomb 
interaction which is close in form with recent work by 
Samaddar and Das Gupta We present and discuss the 
caloric curve, specific heat and phase diagram for undiffer- 
entiated and differentiated protons and neutrons. In this 
framework we work out the cluster content of the system 
and related observables for different values of the tem- 
perature. Most of the simulations are performed in the 
framework of the canonical ensemble. A comparison with 
microcanonical calculations is also presented. In section 
3 we introduce the cellular model and study its thermo- 
dynamic properties which will be confronted with those 
obtained with the lattice model. Section 4 is devoted to 
comments and conclusions. 
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2 Effect of the Coulomb interaction on finite 
systems described by Ising-type lattice 
models 

2.1 The models and their implementation 

The classical system with A particles located on a cu- 
bic lattice with sites is described by Hamiltonians of 
the form 



<»j> 



1 . 



e^Z,Zi 



(1) 



where Ki [i = 1, ■ • ■ , L'^) labels either a proton (p) or a 
neutron (n). The short range potential VxiKj acts only be- 
tween pairs of nearest neighbour particles < i,j >, when 
the sites i and j are both occupied (n^ = Uj = 1), unoccu- 
pied sites correspond to nk — 0. The interaction strengths 
are chosen in two different ways, either Vnn — Vpp — Vnp — 
e = 5 MeV (undifferentiating potential, case U in the se- 
quel) or Vnn = Vpp = 0, Vnp = tnp = 5.33 MeV (dif- 
ferentiating potential, case D in the sequel). The latter, 
a priori more realistic choice has been used by different 
authors the choice of strength is such that the 

classical ground state energy of finite nuclear systems is 
reproduced [p^ . 

The Coulomb term in (|^) is such that = or 1 
if particle j is a neutron or a proton respectively. The 
distance between particles i and j, rij, is determined on 
a lattice in which the distance between sites is fixed to 
d= 1.8 fm. 

The canonical partition function reads 



Z(/3) = 



E 



■S- 



(2) 



where Z is the total number of charges in the system. In 
the calculations Z/A was fixed at the value of 0.4. 

The thermodynamic properties and the space occupa- 
tion by particles and bound clusters are obtained from re- 
alisations generated by means of Metropolis Monte Carlo 
simulations in which particles are moved on the lattice. 
Technical details about the algorithm used in the frame- 
work of the canonical ensemble have been given elsewhere 
]lO| . We selected 5 x 10'' realisations corresponding each 
to 10 X Metropolis steps. 

We have also performed calculations in the framework 
of the microcanonical ensemble. A given realisation (r) 
with fixed energy E has a weight 



Weir) ^{E- U{r)) — -^ ■ 0{E ~ U{r)) 



(3) 



where U{r) is the potential energy, O is the step function, 
D the space dimension, here D = 3. Detailed balance fixes 
the acceptance rate from a realisation (r) to a realisation 
(r') to 



Wr—>r' = niin 



1, 



Wsir) 



(4) 



The simulations were done with open boundary con- 
ditions on the edges of the lattice both in the case of the 
canonical and microcanonical ensembles. 



2.2 Thermodynamics and cluster size distributions 

The quantities which we consider here are the caloric 
curve and the specific heat for fixed volume which in the 
canonical ensemble reads 



d<U> I , ^ . 



The brackets stand for an average over an ensemble of 
systems. In the microcanonical ensemble the temperature 
is related to the kinetic energy K of the system by 

T = 2< K > /AD . 

Then the specific heat can be cast in the simple form 



DA 



DA (DA A / 1 \' 



DA 



In the absence of the Coulomb interaction clusters 
of bound particles are determined in the same way as 
in ref. H by application of the Coniglio-Klein prescrip- 
tion p7[| which fixes the condition under which a particle 
located on a site which is topologically connected to a 
cluster does effectively belong to this bound cluster. The 
test between nearest neighbour pairs of particles is made 
by means of the probability p = 1 — exp(— e/2T) in case 
U, with e replaced by gpp, e„„, e„p in case D, depending on 
the nature of the particles which are involved. The prob- 
ability p is then compared to a random number ^ G [0, 1]. 
If ^ ^ P the particles are considered as being bound, if 
f > p they are not. 

The Coulomb interaction is taken into account in the 
case U and for a pair of neighbouring charged particles 
by replacing e by e — e^/d where d is the (fixed) dis- 
tance between the particles. In the case D p is not mod- 
ified since neighbouring charged particles are not bound 
Epp = 0. This procedure deserves some comments. The 
identification of fragments by means of the Coniglio-Klein 
prescription is conceptually valid if the interaction is short 
ranged, i.e. concerns only nearest neighbours. This is not 
the case with the Coulomb interaction which acts over 
a large range. We restricted the tests with Coulomb for 
the connection between particles to nearest neighbours. 
There is no real justification for this, since it is not possi- 
ble to invoke screening effects, the classical Debye-Hiickel 
screening length being much larger than the size of the 
system. Hence, as it stands, we have no real means to 
estimate the validity of our procedure, except for quan- 
titative arguments, i.e. the effect is certainly weak since 
e^/c? ~ 0.8 MeV should be compared to e ~ 5 MeV and 
further the Coulomb interaction acts only on protons in a 
system where N/Z =1.5 {N = number of neutrons). 

The cluster size distribution can be determined for any 
fixed density p ~ A/L^ and temperature T. From its 
knowledge one can work out different observables. Here 
we restrict ourselves to the behaviour of the largest clus- 
ter Arnax as a function of log 5*2, with So — m2/mi and 
nik the k*'^ moment of the distribution |l§,[l9]. 
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2.3 Effect of the Coulomb interaction on the 
thermodynamic properties of the system 

The effect of the Coulomb interaction has been studied 
in both cases U and D with the two-body strengths given 
in 2.1. Fig. |l| shows the caloric curve and the specific heat 
associated to a system with L = 10 and p = 0.3 for case 
U. Calculations have been performed in both the canon- 
ical and microcanonical framework. As it can be seen in 
Fig. 0, on the specific heat, the results are undistinguish- 
able. The point of interest is the fact that the Coulomb 
interaction induces a sizable shift in the energy for fixed 
temperature on the caloric curve which increases approx- 
imately parallel to each other and a reduction of about 
1.5 MeV in the temperature corresponding to the maxi- 
mum of Cy- The observed reduction of the temperature 
can be understood by means of the following arguments. 
The Coulomb interaction takes a part of the total available 
energy and hence the thermal energy is decreased leading 
to a decrease of the temperature. With changing excita- 
tion energy the Coulomb energy does not change much, 
hence the energy differences between the cases without 
and with Coulomb interaction stay approximately con- 
stant. This behaviour is general, whatever the size of the 
system and its density. The expected drop in the tem- 
perature in the presence of Coulomb reflects also in the 
(p, T) phase diagram shown in Fig. ||, in which the coexis- 
tence line has been fixed by following the maximum of Cy. 
One may mention the slight dissymmetry with respect to 
p = 0.5. For larger p the Coulomb effect is larger than 
for small densities and consequently produces a stronger 
shift in the temperature which defines the border of the 
coexistence region. This is due to the fact that higher p 
corresponds to stronger packing, hence charges interact 
more effectively. In case D, when the particles are differ- 
entiated, the results are qualitatively similar for both the 
caloric curve and the specific heat with a similar shift of 
1.5 MeV in the temperature, see Fig. ^ In this case, the 
maximum of Cy is much more damped when Coulomb is 
present than in case U, and the transition region is there- 
fore less precisely defined. 

One may notice that in both cases, with and without 
Coulomb interaction, the caloric curve is unable to give 
an indication about the order of the transition. In former 
studies [|lO| it has been shown through the attempt to 
determine critical exponents that the transition could be 
continuous, but it came out that the result was not conclu- 
sive because of the small size of the system used in order 
to apply finite size scaling arguments. Microcanonical cal- 
culations on small systems ||2^ also seemed to indicate 
the continuous nature of the transition, no backbending 
in the caloric curve could be observed. This point was 
later discussed in rcf. [| l 2 | and recently it was shown by 
Pleimling and Hiiller pl| that the backbending can be ob- 
served if the system gets large enough, which shows that 
the transition is in fact first order when one goes to the 
thermodynamic limit. 



2.4 The effect of the Coulomb interaction on the 
cluster size distributions and related observables 

The behaviour of the mass distribution of clusters is 
shown in Fig. ^ for case D, temperatures lie in the range 
3 — 5 MeV. One observes again a shift in the behaviour 
of the system when the Coulomb interaction is switched 
on. For the same temperature, the distribution without 
Coulomb contains more larger clusters. This is under- 
standable since the long range interaction is repulsive and 
tends to split the bond fragments. The whole picture is 
also consistent with the thermodynamic behaviour of the 
system. This is seen in Fig. ^ on the line which separates 
the heavy cluster behaviour below the line from the light 
cluster behaviour above it and stands as the finite size 
remnant of a continuous set of transition points in the 
percolation framework |2^. Indeed, the line lies at lower 
temperatures when the Coulomb interaction is present. 
The Coniglio-Klein procedure is the correct prescription 
in order to identify bound clusters of interacting parti- 
cles with kinetic energy. It ensures that a given particle is 
bound to a cluster or not by means of separation energy 
arguments Q . One would expect that this separation line 
ends at p = 0.5 on the coexistence line. This is not the 
case here. It may be due to the criterion (the maximum 
of Cv) which is used here in order to fix this coexistence 
line. It clearly corresponds to a finite size effect. Indeed 
we checked that the distance between the separation line 
and the p = 0.5 point on the coexistence line gets smaller 
and smaller when the size of the system increases. 

Fig. 1^ shows the correlation between Amax and ^2 in 
case D. Events which come at large Amax correspond to 
the presence of heavy clusters, those with small Amax to 
the presence of many small clusters. The intermediate re- 
gion shows events which correspond to the transition re- 
gion in the percolation framework. As it can be seen, this 
region appears again for lower temperatures in the case 
where the Coulomb interaction is acting. 

A further observable of physical interest which can 
be experimentally determined is the ratio of the number 
of neutrons to the number of protons, N/Z, which are 
present in the clusters . Fig. ^ shows the evolution of 
this ratio for clusters of mass A = N + Z > with in- 
creasing temperature for case D. One observes a decrease 
of the relative number of neutrons, both in the absence and 
the presence of the Coulomb interaction. This is in agree- 
ment with the results of ref. ||T^ and, for heavy clusters, 
with the experimental results of ref. The decrease 

is apparently slightly more pronounced in the case where 
the Coulomb interaction is active, although, as indicated, 
the fluctuations from event to event increase with increas- 
ing temperature. This result may appear to be somewhat 
paradoxical. It may be due to the fact that for fixed tem- 
perature the cluster size distribution is different in both 
cases as we already saw above. 

Finally, the evolution of the N/Z ratios with the mass 
is shown in Figs. |^ for case D. In Fig. ^ the results cor- 
respond to the case where the Coulomb interaction is in- 
cluded. The ratio drops from N/Z ^ 1.6 for A = 3 to 1.2 
for the largest masses. A similar behaviour is observed in 
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Fig. ^ for a higher temperature and in the absence of 
the Coulomb interaction. Hence lighter species are more 
neutron rich than heavier ones. This fact which has al- 
ready been observed in the calculations of ref. @ in the 
absence of the Coulomb interaction remains valid when 
Coulomb is present and seems to be in agreement with 
the experimental findings [p5|,p6[. 



3 Cellular model approach to the description 
of fragmenting nuclei 

A system of A particles is contained in a cube of 
volume V in 31? space and divided into cells of volume 
d^, V = L^(P. Cells are either empty or occupied by 
at most one particle characterized by its random posi- 
tion Ti (j = 1, • • • , yl) in the cell and random momen- 
tum Pi 10]. The particles are classical. Neighbouring par- 
ticles interact through a short range two-body potential 
VQ{rij) which is repulsive at short distance and attractive 
for Tij > 1.55 fm In the calculations d = 1.8 fm. 
The Coulomb interaction is not taken into account. The 
Hamiltonian is written as 



P. 



=1 



(5) 



with rii = 0(1) if the cell i is empty (occupied). 

Similarly to the lattice case we have worked out the 
thermodynamic properties of such a system in the frame- 
work of both the canonical and microcanonical ensemble. 
The generation of configurations by means of a Metropo- 
lis Monte Carlo procedure is performed in the following 
way. Starting from an initial configuration in which the 
particles are disposed randomly into the available cells 
one performs either the move of a particle into an empty 
cell, or changes the position of a particle in its cell. Each 
operation is performed with a probability 1/2, the appli- 
cation of the Metropolis Monte Carlo algorithm leads to 
a minimum in energy in the framework of the canonical 
ensemble. In the case of the microcanonical ensemble the 

(i- 



Eq. 



test is the one described in section p.l| , ^ 

Typical results can be seen on Fig. || which shows the 
caloric curve and the specific heat. Both ensembles lead to 
a very close behaviour of these quantities. This behaviour 
is qualitatively similar to those obtained with lattice mod- 
els. The same is true for the corresponding phase diagram 
shown in Fig. ^ The thermodynamic phase coexistence 
line shows two slight maxima symmetrical with respect to 
the critical point at the density p — 0.5. This behaviour 
may be due to the way in which we define the coexistence 
line (maximum of the specific heat) and (or) related to fi- 
nite size effects, see ref. [|o). The present calculations have 
been performed for a fixed volume with 8^ sites and the 
possible transition looks continuous in this case. The small 
depression for p = 0.5 reminds the same effect which was 
found in the IMFM calculations of ref. 0. The dotted 
line shows the separation line of heavy and light clusters 
as in the case of the lattice models. Its behaviour is very 



similar to the one observed in Fig. The fact that its 
lower end does not coincide with p = 0.5 m ay have the 
same reasons than those presented in section 2.4. 



4 Summary and conclusions 

In the present work we aimed to present and discuss 
two points related to the description of nuclear fragmen- 
tation by means of microscopic classical models. 

We first looked for the effects induced by the long range 
Coulomb interaction which acts between charged particles 
in the presence of a short range attractive potential mim- 
icking the nuclear interaction. We have shown on a lattice 
model that the Coulomb interaction does not induce any 
qualitative change in the thermodynamics and cluster size 
distribution in the different phases in which the system 
exists. The essential quantitative effect is a global system- 
atic and sizable shift of the temperature of the system by 
about 1.5 MeV in the models which have been worked out 
both in the thermodynamic and cluster transitions. The 
same type of effect has already been mentioned in former 
studies We have also investigated the behaviour of 
the cluster content, both in the absence and presence of 
the Coulomb interaction. 

In a second part we introduced a cellular model aimed 
to describe a disordered system of particles in thermo- 
dynamic equilibrium in the framework of the canonical 
and microcanonical ensembles. The results show that such 
models, though a priori more realistic than lattice mod- 
els, are descriptions which lead qualitatively to the same 
properties, at least when finite systems are concerned. 
This leads to the conclusion that lattice models which are 
simple to handle should be good enough for a schematic 
description of nuclear fragmentation processes if the gen- 
erated systems are in thermodynamic equilibrium. There 
remains however the problem concerning the importance 
of quantum effects. These may be the weaker the higher 
the temperature, but this point has to be confirmed. 
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Fig. 1. Caloric curves and specific heat for a lattice model with 
A = 300 undifferentiated particles (protons and neutrons) in 
a volume = 10^. Full and open lines correspond respec- 
tively to results including or not the Coulomb interaction. The 
specific heat calculated with the Coulomb interaction in the 
framework of the microcanonical ensemble is represented by a 
dotted fine. 
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Fig. 2. Phase diagrams (p, T) with Coulomb interaction 
(left) and without Coulomb interaction (right) for undifferenti- 
ated particles. The full line indicates the phase separation, the 
dashed line the separation line between systems with heavy 
clusters (below) and light clusters and particles (above). See 
text. 
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Fig. 3. Same as Fig. |ljfor differentiated particles (protons and 
neutrons). Calculations are done in the canonical ensemble. 
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Fig. 4. Cluster size multiplicities for different temperatures. Fig. 5. Size of the heaviest cluster Amax vs. S2 = m2/m\ rep- 
Left column: with Coulomb interaction. Right column: without resented for the different Monte Carlo events. m2 and mi are 
Coulomb interaction. See text. the second and first moments of the cluster size distribution. 
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Fig. 6. Average isotopic ratios N/Z and fluctuations widths Fig. 7. Isotopic ratios N/Z in the clusters as a function of the 

in the clusters as a function of the temperature in the presence cluster size. Upper part (a): with Coulomb interaction. Lower 

and absence of the Coulomb interaction. Volume = 10^, part (b): without Coulomb interaction, 
density p = 0.3. See comments in the text. 
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Fig. 9. Phase diagrams (p, T) in the framework of the ceUuIar 
model for a system with A — 256 particles in a volume = 8"^ . 
The dashed line has the same meaning as in Fig. ^ 
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Fig. 8. Canonical and microcanonical caloric curves and spe- 
cific heat in the framework of the cellular model for a system 



of A = 
text. 



256 particles in a volume L = 8 . See comments in the 



